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We present the Large Hadron Collider (LHC) discovery potential in the Z' and 
heavy neutrino sectors of a U{1)b-l enlarged Standard Model also encompassing 
three heavy Majorana neutrinos. This model exhibits novel signatures at the LHC, 
the most interesting arising from a Z' decay chain involving heavy neutrinos, even- 
tually decaying into leptons and jets. In particular, this signature allows one to 
measure the Z' and heavy neutrino masses involved. In addition, over a large region 
of parameter space, the heavy neutrinos are rather long-lived particles producing 
distinctive displaced vertices that can be seen in the detectors. Lastly, the simulta- 
neous measurement of both the heavy neutrino mass and decay length enables an 
estimate of the absolute mass of the parent light neutrino. 
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I. INTRODUCTION 

The B — L (baryon number minus lepton number) symmetry plays an important role in 
various physics scenarios beyond the Standard Model (SM). Firstly, the gauged U{1)b-l 
symmetry group is contained in a Grand Unified Theory (GUT) described by a 5*0(10) 
group [1] . Secondly, the scale of the B — L symmetry breaking is related to the mass scale 
of the heavy right-handed Majorana neutrino mass terms providing the well-known see-saw 



mechanism {2] of light neutrino mass generation. Thirdly, the B — L symmetry and the scale 
of its breaking are tightly connected to the baryogenesis mechanism through leptogenesis 
via sphaleron interactions preserving B — L. 

In the present paper we study the minimal B—L low-energy extension of the SM consisting 
of a further U{1)b~l gauge group, three right-handed neutrinos and an additional Higgs 
boson generated through the U{1)b-l symmetry breaking. It is important to note that in 
this model the B — L breaking can take place at the Electro- Weak (EW) or TeV scale, i.e., 
a value far below that of any GUT scale. This B — L scenario therefore has potentially 
interesting signatures at hadron colliders, particularly the LHC. New particle states such as 
Z', Higgses and neutrinos, all naturally have masses at the EW or TeV scale. The breaking 
of the B — L symmetry at the EW or TeV scale can be viewed as a remnant of a grand 
unified gauge symmetry, such as 5*0(10). Furthermore, with respect to baryogenesis, since 
B + L IB, violated by sphaleron interactions, this implies that baryogenesis or leptogenesis 
cannot occur above the scale of B — L breaking. A scenario with B — L breaking at the EW 
(or TeV) scale therefore implies EW (or TeV) scale baryogenesis j^. 

The particular subject of the present paper is the first detailed study of the collider 
phenomenology of the gauge and fermionic sectors of the minimal B — L extension of the 
SM, where the additional U{1)b-l gauge group is indeed associated to the B — L number 
|l, ^,[6]. The analysis of the scalar sector will appear in a future paper f?]. The new results 
on B — L phenomenology at the LHC include observable signals from a Z'-boson as well as 
heavy neutrinos with a mass of up to several hundred GeV^. A very interesting feature of 
such a B — L model is possibly relatively long lifetimes of the heavy neutrinos which can 
directly be measured. In turn, such a measurement could be a key to sheding light on the 



^ We deliberately assume that the B~ L symmetry breaking scale generating Mz' is somewhat higher than 
the scale of the heavier right-handed neutrino mass, thereby enabling Z' to heavy neutrino decays. 
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mass spectra of the light neutrinos. 

This work is organised as follows. Sect. HTl reviews the model under study and its imple- 
mentation in the CalcHEP package 8|], together with an overview of its parameter space. 
Sects. IIII Al and IIIIBI study the decay properties - specifically, width and Branching Ratios 
(BRs) - of the new spin-1 and spin-1/2 particles of the B — L model. Sect. IIII CI discusses 
their experimental signatures and production plus decay cross sections and also contains 
a numerical analysis for two particular benchmark points in the B — L parameter space, 
followed by a study of the expected background. The conclusions are in Sec. HVl 



II. THE B - L MODEL AND ITS IMPLEMENTATION INTO CALCHEP 

A. The model 

The model under study is the so-called "pure" or "minimal" B — L model (see 
for conventions and references) since it has vanishing mixing between the two U{1)y 
and U{1)b-l groups. In the rest of this paper we refer to this model simply as the 
"5 — L model". In this model the classical gauge invariant Lagrangian, obeying the 
SU{?>)c X SU{2)l X U{1)y X U{1)b-l gauge symmetry, can be decomposed as: 

= ^YM + + -^f + -^Y ■ (1) 

The non-Abelian field strengths in J^ym are the same as in the SM whereas the Abelian 
ones can be written as follows: 

where 

Ffj.u = d,_tBy — dyB^ , (3) 
= d,Bl - d^B'^ . (4) 

In this field basis, the covariant derivative is: 

D^ = d^ + igsT'^G^^ + igT'^W;' + igiYB, + ig[YB-LB'^ . (5) 
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The fermionic Lagrangian (where k is the generation index) is given by 

3 

k=l 

+ilkL'lt,D^'lkL + ie^lf^D'^ekR + iWRjf.D'^UkR^ , (6) 

with the respective fermion charges given in Tab. [H The B — L charge assignments of new 
fields as well as the introduction of new scalar Higgs (x) and fermionic right-handed heavy 
neutrinos {vr) fields are designed to eliminate the triangle B — L gauge anomalies. (Tab. HTl 
shows the scalar content and charges of our B — L model.) Therefore, the B — L gauge 
extension of the SM group broken at the EW scale does necessarily require at least one 
new scalar field and three new fermionic fields which are charged with respect to the B — L 
group. 

The scalar Lagrangian is: 

if, = {D^^H)^ D^H + {D^x)^ D^x - V{H, x) , (7) 

with the scalar potential given by 

V{H, x) = m'H^H + /x^ I X I' +Xi{H^H)' + X2 \ X \' +\,H^ H \x\\ (8) 

where H and x are the complex scalar Higgs doublet and singlet fields, respectively. 
Finally, the Yukawa interactions are: 

■^Y = -y%qjLdkRH - yJkOjL^kRH - yji^ljLCkRH 

-y'^jkhLT^kRH - vfkMji^kRX + h.c. , (9) 

where H = ia^H* and i,j,k take the values 1 to 3, where the last term is the Majorana 
contribution and the others the usual Dirac ones. 



B. Model implementation into CalcHEP 



We make use of the CalcHEP package [8| to study the collider phenomenology of the B — L 
model. For the derivation of the Feynman rules (see the Appendix for those pertaining to 
the heavy neutrino interactions) and for the straightforward implementation of the model in 



the CalcHEP package, we have used the LanHEP module 



3 . The availability of the model 
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TABLE I: Fermion content and charges for the B — L model. 
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TABLE IL Scalar content and charges for the B — L model. 

implementation into CalcHEP in both the unitary and t'Hooft-Feynman gauges allowed us 
to perform powerful cross-checks to test the consistency of the model itself. 

The implementation of the gauge sector is quite straightforward. Since there is no mixing 
between the (SM) Z and Zb-l bosons (hereafter, we will refer to the Zb~l boson as a Z') one 
just needs to define a new heavy neutral gauge boson together with the covariant derivative 
given by eq. ([5]) and the charge assignments in Tabs. [114111 For the scalar sector, we need to 
implement the mixing between mass and gauge eigenstates of the two Higgs bosons. 

The implementation of the neutrino sector is somewhat more complicated. Majorana- 
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like Yukawa terms are present in eq. ([9]) for the right-handed neutrinos, therefore one must 
implement this sector such that the gauge invariance of the model is explicity preserved. 
This can be done as follows. As a first step we rewrite Dirac neutrino fields in terms of 
Majorana ones using the following general substitution: 

D 1-75 1 + 75 
^ = -^—^L + -^—^R ' (10) 

where is a Dirac field and i^l{r) are its left (right) Majorana components. If we perform 
the substitution of eq. (|TOl) in the neutrino sector of the SM, we will have an equivalent theory 
formulated in terms of Majorana neutrinos consistent with all experimental constraints. 
The second step is to diagonalise the neutrino mass matrix from eq. ([9]): 

tud , 

, (11) 



niD M 



where 



M = V2y^x, (12) 

where x is the Vacuum Expectation Value (VEV) of the x field. This matrix can be diago- 
nalised by a rotation about an angle a^, such that: 

tan 2a,, = — —. (13) 

For simplicity we neglect the inter-generational mixing so that neutrinos of each gen- 
eration can be diagonalised independently. We also require that the neutrinos be mass 
degenerate. Thus, z/^^r can be written as the following linear combination of Majorana mass 
eigenstates ui^h '■ 

(14) 

The last subtle point is the way the Lagrangian has to be written, in particular the 
Majorana-like Yukawa terms for the right-handed neutrinos (the last term in eq. ([9])). In 
order to explicitly preserve gauge invariance, this term has to be written, in two-component 
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notation, as: 



- /'^^^^^X + h.c. , (15) 
where v is the Dirac field of eq. ffTOl) . whose Majorana components i'l,r mix as in eq. ([1 
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C. Parameter space 

In this section we define the independent parameters of the B — L model and their vahd 
range^. The set of parameters relevant to our study is the following one. 

• the new U{1)b-l gauge coupling. Here, the absence of a Landau pole up to the 
GUT scale Qgut = 10^^ GeV implies g'l < 0.5 from a Renormalisation Group Equation 



(RGE) analysis %m\- 



Mz'i the new gauge boson mass. An indirect constraint on Mz' comes from analyses at 
;he Large Electron-Positron (LEP) collider of Fermi effective four-fermions interactions 



ll|: 



^>6TeV. (16) 
9i 

As we demonstrate below, this constraint will provide an upper bound to the Z' 
production cross sections at the LHC. 

• M,^^, the heavy neutrino masses. We take them to be degenerate and relatively light, 
varying in the range 50 GeV < M^^ < 500 GeV. 

• mj,;, the SM (or light) neutrino masses. We use the cosmological upper bound 

mu^ < 1 eV. As we will see in Sect. IIII A[ detectable displaced vertices may 
occur for rriyi < 10^^ eV. 

For this analysis of the collider phenomenology of our B — L model we have chosen 
= 1.5 TeV and g'^ = 0.2 as a representative point satisfying present experimental 
constraints as well as two values for heavy neutrino masses: M,y^ = 200 and 500 GeV. 
Finally, we have fixed the light neutrino mass to be 171,^^ = 10~^ eV. For illustrative purposes 
we take all neutrino masses, both light and heavy, to be degenerate. 



^ Since the scalar sector is not within the scope of the present study we do not discuss the corresponding 
parameters, as one can choose settings in parameter space such that the scalars are entirely decoupled 
from the remaining particles. Here, we achieve this by requiring Ai = 3, A2 = 0.08 and A3 = 0.01, so 
that the scalars masses are mh-^ « 600 GeV and to/i^ w 1.5 TeV, corresponding to our default benchmark 
parameters in the Z' sector, g[ — 0.2, Mz' — 1.5 TeV. The only exceptions to this will be made in 
Sects. IIIA and IIIIBi where we will temporarily adopt other settings, to describe the complete decay 
pattern of the Z' and heavy neutrinos. 
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III. PHENOMENOLOGY OF THE B~L MODEL 
A. Heavy neutrino properties 

1. Heavy neutrino decays 

After the diagonalisation of the neutrino mass matrix reahsing the see-saw mechanism, 
we obtain three very hght neutrinos (ui), which are the SM-hke neutrinos, and three heavy 
neutrinos (i^h)- The latter have an extremely small mixing with the ui^s thereby providing 
very small but non-vanishing couplings to gauge and Higgs bosons (see the Appendix for 
the Feynman rules involving heavy neutrino interactions) which in turn enable the following 

decays: Uh — > l^W"^, —>■ uiZ, —>■ i^ihi, Vh ^ih2 as well as — > i^iZ' when these 
decay channels are kinematically allowed. Fig. [1] presents the corresponding BRs versus 
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FIG. 1: Heavy neutrino branching ratios versus its mass for the fixed Mz' = 1-5 TeV, m^j^ = 150 
GeV and m/,^ = 450 GeV corresponding to Ai = 0.19, A2 = 0.017, A3 = 0.01 and g[ = 0.2. Here, 
W means the sum over and W~ . 

the heavy neutrino mass for the values of the other relevant parameters given in the 
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caption. One can see that the BRiyh l^W^) is dominant and reaches the 2/3 level in 
the My^ ^ My/,Mz limit, while BR{vh i^iZ) and BR{uh —>■ vihi) both reach the 1/6 
level in this regime. In contrast, the Vh —>■ ^«^2 as well as Uh viZ' decay channels are well 
below the percent level and are negligible for our study. In this paper we will eventually 
assume that the heavy neutrino masses are smaller than both Higgs boson masses. Under 
this assumption — > vihi [i = 1,2) is not kimematically possible and BR{vh — > I'lZ) 
reaches the 1/3 level in the m^^ ^ My/^ Mz limit. 



2. Lifetime of the heavy neutrinos 

The heavy neutrino couplings to the weak gauge bosons are proportional to the ratio of 
light and heavy neutrino masses (see the Appendix), which is extremely small. Therefore the 
decay width of the heavy neutrino is correspondingly small and its lifetime large. The heavy 
neutrino can therefore be a long lived particle and, over a large portion of parameter space, 
its lifetime can be comparable to or exceed that of the 6-quark. (In fact, for my^ = 10~^ eV 
and = 200 GeV they are equal.) In Fig. [2] we present the heavy neutrino lifetime (top) 
in pico-seconds and the proper decay length (or mean path) (bottom) in micro-meters as a 
function of the light neutrino mass. The proper decay length is defined as ctq, where tq is 
the lifetime of the heavy neutrino. The purple band presents the proper decay length of the 
6-quark while the blue band indicates the range of a typical micro-vertex detector. The red 
band shows the region of light neutrino masses excluded by direct measurements of neutrino 



oscillations |12)], by taking the lighter neutrino to be massless (so that the other neutrinos 
cannot populate this region). One should also note that the lifetime and the proper decay 
length of the heavy neutrinos in the laboratory frame will actually be equal to those given in 
Fig. [2] times the Lorentz factor equal to p^j^/M^^^ defined by the ratio Mz> jMy^ which can be 
as large as about a factor of 10. We can then see that there exists a region where the heavy 
neutrino lifetime is of the same order as that of the 6-quark (shown as a purple band). The 
mean path and the respective lifetime of heavy neutrinos can therefore be measured from a 
displaced vertex inside the detector. The heavy neutrino can however be distinguished from 
a 6-hadron through the observation of vertices consisting of only two isolated leptons. (A 
SM S-meson decay while possible would have a very small BR, ~ 10~^ at the most.) 

An experimentally resolvable non-zero lifetime along with a mass determination for the 
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FIG. 2: Heavy neutrino life-time (top) and proper decay length (or mean path) ctq (bottom) as 
a function of the light neutrino mass. The purple band presents the proper decay length of the 
6-quark while the blue band indicates the range of a typical micro-vertex detector. The red band 
shows the region excluded by neutrino oscillation direct measurements. 

heavy neutrino also enables a determination of the light neutrino mass. The lifetime mea- 
surement allows the small heavy-light neutrino mixing to be determined and as one can see 
from eqs. f|TT]) - f|T^ this, along with the heavy neutrino mass, gives the light neutrino mass. 
Considering only one generation for simplicity, this is expanded upon below. 

Mass eigenstates are related to gauge eigenstates by eq. ffT^ . hence the eigenvalues are 
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given by solving the equation: 



M, 




rriD 



M 



which yields 



m,. 



sin 2au ttid + sin^ M . 



= — sin 2a,,m£) + cos^ M . 



(17) 
(18) 



with ftj, given by eq. f|T3l) . We have then three parameters {rriD, M and a^) and a constraint 
(given by eq. (fT3l) ). that can be used to eliminate one parameter from the above equations. 

The Feynman rules given in the Appendix demonstrate that heavy neutrino interactions 
are determined by the mixing angle only, as is the total width (and therefore the mean 
decay length). Hence, it is convenient to keep M^^ and a^, as independent model parameters 
eliminating from eq. f|T3l) . 

rriD = moioiu, M^J . (19) 
By measuring the heavy neutrino mass we can also invert eq. (fTSll 



(20) 



to finally get a fully known expression for the SM light neutrino mass as a function of 
our input parameters M^^^ and a^, that we can measure independently, by inserting eqs. 
([T9])-(l20l) into eq. (HT]), 

m^X^D,M) = m^^{a^,M^J . (21) 

This simple picture shows that within the B — L model we have an indirect way of accessing 
the SM light neutrino mass by measuring the mass of the heavy neutrino and the kinematic 
features of its displaced vertex. If the whole structure of mixing is taken into account, 
including inter-generational mixing in the heavy neutrino sector, the task of determining the 
light neutrino mass this way would become more complicated but the qualitative features 
and the overall strategy would remain the same thereby providing one with a unique link 
between very large and very small mass objects. 
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B. Z' decay properties 

As discussed earlier, the extra U{X)b-l gauge group provides an additional neutral gauge 
boson, Z', with no mixing with the SM Z-boson. Therefore our Z' boson decays only to 
fermions at tree-level and its width is given by the following expression: 



2 

1 + 2^ 



4mf 

Mi:- (^^' 



where m/ is the mass and C/ the number of colours of the fermion type / and 
{B — L) X g[ is the vector coupling (see Tab. [T]). 



In Figs. 3(a) and 3(b) we present the total decay width of the Z' as a function of Mz' and 



g'l, respectively (with the other parameters held fixed to three different values), assuming 



that the partial decay width into heavy neutrinos vanishes. Also, Fig. 3(c) presents the 
relative variation of the total width as a function of the mass for three different values of 
Mz' and with g[ = 0.5. 

From the first two plots we see that the total width of a Z' gauge boson varies from a few 
to hundreds of GeV over a mass range of 0.5 TeV < M^/ < 5 TeV, depending on the value 
of g'l, while from the third plot one can gather the importance of taking into consideration 
the heavy neutrinos, since their relative contribution to the total width can be as large as 
25% (whenever this channel is open). One should also note that possible Z' decays into 
one light and one heavy neutrino are highly suppressed by the corresponding (heavy-light) 
neutrino mixing and thus they can safely be neglected. 

The possibility of decays of the Z' gauge boson into pairs of heavy neutrinos is one of 
the most significant results of this work since, in addition to the clean SM-like di-lepton 
signature, it provides multi-lepton signatures where backgrounds can strongly be supressed. 
In order to address this quantitatively, we first determine the relevant BRs. Clearly, these 
depend strongly on the heavy neutrino mass and Fig. H] shows how they change with fixed 
(although arbitrary) values of M^^, for the following three cases: a heavy neutrino (i) much 
lighter than, (ii) lighter than and (iii) comparable in mass to the Z', in the range 0.5 TeV 
< Mz' < 5 TeV, before summing over generations. 

A feature of the current B — L model illustrated in the previous figures is that the Z' 
predominantly couples to leptons. In fact, after summing over the generations, k = 1...3, we 
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FIG. 3: Z' total width as a function of: 3(a) Mz' (for fixed values of g^), |3(b)"] g^ (for fixed values 
of Mz') and|3(c)|M^^ (for fixed values of Mz' and g[ = 0.5). 
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FIG. 4: Z' BRs as a function of Mz' for several heavy neutrino masses: M^^ = 50, 200 and 1000 
GeV, from left to right, respectively. A summation over all lepton/neutrino flavours is implied 
throughout whereas in the case of quarks we distinquish between light flavours {q = d,u,s,c,b) 
and the top quark. 

roughly get for leptons and quarks: 

^ BR [Z' Iklk + Vki^) ~ I , ^BR{Z' ^ qkW^) ~ - . 

k k 

Not surprisingly then, for a relatively light (with respect to the Z' gauge boson) heavy 
neutrino, the Z' BR into pairs of such particles is relatively high: ~ 18% (at most, again, 
after summing over the generations). 

Combining this last result together with those of Sect. IIII Bl we can discuss an interesting 
feature of the B—L model, namely, the multi-lepton signatures (meaning two or more leptons 
being involved). A single heavy neutrino decay will produce a signature of 0, 1 or 2 charged 
leptons, depending on whether the heavy neutrino decays via a charged or neutral current 
and on the subsequent decays of the SM and Z gauge bosons. We can have both chains 

Uh^l^W^ ^l^ + I (23) 

I hadrons 

and 

Uh^uiZ^ui+l . (24) 

I z// z/^/hadrons 

The pattern in ( l23l) provides 1 or 2 charged leptons whilst that in ( IMl) zero or 2, so that 
multi-lepton signatures may arise when the Z' gauge boson decays into a pair of heavy 
neutrinos, producing up to four charged leptons in the final state. Fig. [5] shows the BRs of a 
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FIG. 5: Z' BRs, as a function of M^^, into: 2 leptons (both e and /j,, top-left); 3 and 4 leptons + 
X (both e and fi, top-right); 2 leptons -|- X jets (both e and /i, bottom-left); zoom of the previous 
plot with same legend (bottom-right). 



Z' decaying into 2 (top-left) and 3 or 4 (top-right) leptons (plus possibly missing transverse 
momentum and/or jets, as appropriate) as a function of M^^, where a lepton can be either 
an electron or a muon and these contributions are summed. While the former are clearly 
dominant the latter are not at all negligible. For Miy± < M^^ < Mz, the — > 
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decay is the only one kinematically possible whereas for M^^^ < My/± the heavy neutrino 
can decay only via an off shell W and is therefore very long lived. For a very massive Z' (2 
TeV < Mz' < 5 TeV) the multi-leptonic BRs are roughly 2.5% in the case of Z' — > 3/ and 
0.5% in the case of Z' 41, for a wide range of heavy neutrino masses. 

Finally, from ( l23l) and (12^ one can see that di-lepton decays are possible, whereas Z' 
decays give rise to 2 leptons plus a large amount of missing transverse momentum and/or 
highly energetic jets: see Fig. [5] (bottom- left). Particularly interesting is the decay into 
2 leptons and 4 jets, since here there is no missing transverse momentum at all and its 
BR is rather large with respect to the other non-SM signatures, as we can see in Fig. [5] 
(bottom- left). 



C. Signal-to-background analysis 

In this section we perform a signal-to-background analysis to check the observability at 
the LHC of some of the signatures discussed that may originate from the present B — L 
model. 

In our model setup, wherein the scalar sector is entirely decoupled, all interesting B — L 
signals come from Z' production, since the Z' is the only new particle whose couplings to 
the SM partons are large. The most efficient hadro-production process involving a Z' boson 
is the Drell-Yan (DY) mode 

qq^Z', (25) 

where q is either a valence-quark or a sea-quark in the proton. At the parton level, the Z' 
production cross section for process (!25ll depends on two main parameters: Mz' and g'^. In 



Fig. 6(a) we present the Z' hadro-production cross section a at the LHC as a function of 



both Mzi and g[, in the ranges 0.5 TeV < Mz' < 5 TeV and 0.1 < g[ < 0.5, respectively. 



while Fig. 6(b) presents the contour levels in the (M^/, g'^) plane for cr = 4 pb, 0.3 pb, 50 



fb and 5 fb. The shaded area in Fig. 6(b) is excluded by eq. f|T6l) . 

We expect that the LHC will discover a Z' boson from our B — L model in the standard 
di-lepton decay channel. In Fig. [7] we therefore show the Z' line-shape, i.e., the differential 
cross section for process (1251) as a function of the invariant mass of its decay products, e.g., 
as obtained from the Z' decay into a pair of muons: 

qq^Z'^fi+fi-, (26) 
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FIG. 6: Z' hadro-production cross section at the LHC over the (M^/, g'^) plane: 6(a 



as a continuous 



function of Mz' and g'^ and |6(b)| in the form of contour lines for four fixed values of production 
rates. The dark shaded area on the right-hand side plot is the region excluded by LEP constraints, 
see eq. (fTUj) . 

for the following values of the input parameters: Mzi = 1.5 TeV, g'l = 0.1 4- 0.5 (in 0.1 
steps) and M^^ = 200 GeV. While the di-lepton mode is a powerful Z' discovery channel, 
its sensitivity to the presence of heavy neutrinos is however only indirect, through the Z' 
width, and in fact very weak, as Tz' varies never more than 20% or so due to the presence 
of the new states (recall Fig. E^c)). 

In contrast, multi-lepton signatures carry the hallmark of the heavy neutrinos as the latter 
enter directly the corresponding decay chains and these are explored here by performing a 
detailed Monte Carlo (MC) analysis at the benchmark point M'^ = 1.5 TeV, g[ = 0.2 and 



M, 



200 GeV. The corresponding total cross section for Z' production and decay into 



heavy neutrinos is 46.7 fb (for CTEQ6L [l3| with = Ml^f. 

Through pairs of heavy neutrinos, other than to fully hadronic decays, which are in- 



When discussing event rates in the following, we will assume an integrated luminosity of ^ = 100 fb ^. 
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FIG. 7: Differential cross sections for qq ^ Z' ^ at the LHC for Mz' = 1500 GeV, 

g[ = 0.1 ^ 0.5 (in 0.1 steps) and M^^ = 200 GeV as a function of M^+^-. 



tractable at the LHC (even in presence of the accompanying missing transverse momentum), 
the Z' can also give rise to 2-, 3- or 4-lepton signatures, for both e and /i in the final state. 
Amongst the latter, we intend to study here the case of 3-lepton decays. The reason is 
twofold. On the one hand, we wish to be able to identify heavy neutrino mediation and the 
presence of only one light neutrino in the 3-lepton mode should enable (transverse) mass 
reconstruction (contrary to the case of the 4-lepton channel, where two light neutrinos are 
involved^ On the other hand, we ought to minimise the impact of large backgrounds, so that 



^ Notice that the 4-lepton final state was discussed in See instead Ref. 15| for a discussion of the 

2-lepton signature. 



19 



we neglect here 2-lepton channels (which could easily by overwhelmed by SM DY and tt 
production). 

When the heavy neutrino decays via the PW^ mode, with a subsequent leptonic decay 
of the W^, the charged pair of leptons can carry an invariant mass equal to or lower than 
the heavy neutrino mass, with the maximum invariant mass configuration occurring when 
the light neutrino is produced at rest, so that the edge in this distribution corresponds to 
the Uh mass. A peak in such a distribution corresponding to the SM-like Z boson, coming 
from the uZ decay mode for the heavy neutrino will also be present inthis distribution. 
The di-lepton invarient mass distribution is given in Fig. [HI the difference in the two 
distributions illustrates the effect of taking tau lepton decays that produce muons or electrons 
into account. 




200 400 600 800 1000 1200 

M^, (L=1,2) (GeV) 

FIG. 8: Invariant mass of the two most energetic leptons in the Z' — > 3/ decay. The dashed line 
refers to data without taking into account the tau lepton. (Here, ^ = 100 fb^^.) 

While the invariant mass distribution can provide some insights into the mass of the 
intermediate objects, this is not the best observable in the case of the 3/-signature, because 
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the final state neutrino escapes detection. A more suitable distribution to look at is the 



transverse mass defined in 



16], i.e., 

2 

/ I T, rn / ■ \ . 7-.0 / • \ . I if> I \ 

mn 



(^^^ M\vis) + P^{vis) + - (p'Tivis) +fr)' , (27) 

where {vis) means the sum over the visible particles. For the final state considered here we 
sum over the 3 leptons and 2 jets. The transverse mass distribution is shown in figure [9] 
where a peak at the Z' mass can be seen. We can also see evidence for the presence of a 
heavy neutrino by just considering the 2 most energetic leptons and the missing transverse 
momentum, since this is the topology relevant to a Vh decay. The results show that this 
transverse mass peak for the heavy neutrino is likely to be the best way to measure its mass. 
Both of these configurations are shown in Fig. O the signature of this model is that both of 
the above peaks occur simultaneously. (The different shape for tau-mediated decays is also 
shown in Fig. [9] .) 




FIG. 9: The transverse mass of the two most energetic leptons (left) and all the visible particles 
(right) in the Z' 31 decay. The dashed line refers to data without taking into account the tau 
lepton. (Here, ^ = 100 fb^^) 

The SM background to the 3-lepton signature was studied using CalcHEP. (For simplicity, 
from now on, we limit ourselves to the case without leptonically decaying r's.) Making the 
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assumption that the Z' peak has already been identified and its mass measured elsewhere 
(as is likely from Z' —>■ 21 decays)^, we show that the peak in the distribution can be 
seen despite the initially large backgrounds. This enables one to extract a value for M^^ . In 
the evaluation of the background we considered three sources (including generation cuts, to 
improve efficiency): 

- WZjj associated production (0-3; = 246.7 fb, / = e,/i; ARjj > 0.5, Pj^^ > ^0 GeV, 

hi,2 \ < 3), 

- ti pair production {a2i = 29.6 pb, I = e,fi (6-quark not decayed); QCD scale = Mt/2 
to emulate the next-to-leading order cross section; no cuts applied), 

- tilu associated production {a^i = 8.6 fb, / = e,/i; QCD scale = > 20 GeV). 

In the case of WZjj associated production, three leptons come from the subsequent 
leptonic decays of the two gauge bosons. This is the main source of background. From ti 
pair production two isolated leptons come from the decay of the produced from top 
decay and one additional, third lepton, could come from semileptonic B-meson decay. This 
lepton though will be not generically isolated, because of the large boost of the 6-quark from 
top-quark decay. We use this fact to suppress ti background. Finally, tilu will produce 
three isolated leptons resulting in a significant background despite the small production 
cross section. 

The first set of cuts we use is designed to impose generic detector angular acceptances, 
lepton and jet transverse momentum minimal thresholds and to provide isolation for leptons 
and jets: 



^ Though notice that in Figs. 13-14 the Z' is well above the background, so that its mass could well be fit 
- independently of process (pS)) - in the present channel. 
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Selection ^1 

Vh,2,3\ < 2-5, 

> 15 GeV, 

> 10 GeV, 

^1.2 > 40 GeV; 

ARij > 0.5 V/ = 1...3, J = 1,2, 

ARi^i, > 0.2 VZ,r = 1...3, 

Ai?,,, > 0.5; (28) 



where 



AR= ^/A^fTA^. 



We evaluate the background at two benchmark points for the signal. For the signal the 



common parameters are: 



Mz' 



1500 GeV, 



0.2, 



m, 



10-' eV 



(29) 



and two heavy neutrino masses are considered: 



M^, = 200 GeV, 



500 GeV. 



(30) 
(31) 



These two benchmark points provide two kinematically very different examples. In the 
first case the heavy neutrinos are much lighter than the Z' producing highly boosted events. 
In the second, their mass is comparable to Mz'/2, hence close to their production threshold, 
resulting in minimal boost. From a merely kinematic point of view, all other cases will be 



somewhere between these two. 



Special care should be devoted to the treatment of the tt background, given its large 
production rates which, however, as previously mentioned, can be eliminated by enforcing 
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a suitable lepton-ject separation. The impact of the first set of cuts on the signals and ti 
background is illustrated in Tab. Illli The ARij requirement is indeed extremely effective 
and reduces this background by a factor of 2 • 10^^. Contrary, the loss of signal due to 
this cut is reasonably small. Also note that the signal events with the smaller boost have a 
higher efficiency for passing the angular isolation cuts. 



Cut 


M^, = 200 GeV 




M^^ = 500 GeV 




tt 




# of events 


Eff. % 


# of events 


Eff. % 


# of events Eff. % 


no cuts 


482.32(10) 


100 


239.30(10) 


100 


1.28 • 10^ 100 


T] cuts 


346.44(7) 


71.8 


170.79(7) 


71.4 


5.1 • 10^ 43.9 


AR+Pt cuts 


68.043(15) 


14.1 


73.668(31) 


30.8 


99.699(3) 0.014 



TABLE III: Efficiencies of the Selection #1 cuts for the two benchmark signals and the ti back- 
ground, for events with three or more leptons and with two or more jets in the final state for 
^ = 100 fb^^. In case of ARjj < 0.5 partons were merged into one 'jet' at the very beginning of 
the selection. 

Figs. [TOl - fTT] show the distributions in Mjj, M^^j^, M^jj and after Selection #1 cuts, 
for the signal with the two heavy neutrino masses, 200 and 500 GeV, that we are considering 
and the and backgrounds. 

In signal events both jets come from the therefore we apply the following constraint: 

Selection #2 

\Mjj - Mw\ < 20 GeV. (32) 

After the application of this cut the other distributions considered are shown in Figs. [T2] - 
[T3] for the 200 and 500 GeV heavy neutrino masses, respectively (here, we now also show 
the difference between the and Mijj distributions). From these plots it is clear that 
transverse mass M|^ - • provides good discrimination between signal and background. The 
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FIG. 10: Signal = 200 GeV) and background distributions after the Selection #1 cuts. (Here, 
^ = 100 fb-i.) 

following cut is then used to further suppress the background. 

Selection ^3 

\M^i2j -Mz'\< 250 GeV. (33) 
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FIG. 11: Signal = 500 GeV) and background distributions after the Selection #1 cuts. (Here, 
^ = 100 fb-i.) 

After this set of cuts we end up with a very clean signal for both a 200 and 500 GeV Uh 
mass in the di-lepton transverse mass distribution, in fact practically free from background, 
as shown in Fig. [TH Notice that this variable was formed from the two closest (in ARu) 
leptons since they are likely to originate from the same boosted (unlike Fig. M ^ where 
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FIG. 12: Signal (M^^ = 200 GeV) and background distributions after the Selection ^1 and ^2 
cuts. (Here, ^ = 100 fb^i.) 

the two most energetic leptons were used.) 

In order to establish the signal we finally select events around the visible peak, by 
requiring: 
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FIG. 13: Signal {M^^^ = 500 GeV) and background distributions after Selection ^1 and ^2 cuts. 
(Here, ^ = 100 fb^^.) 

Selection ^4 

0< M^i< 250 GeV or 400 GeV < M^i < 550 GeV, (34) 
depending on the benchmark signal under consideration. The efficiencies of the Selection 
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FIG. 14: Signal {M^^ = 200 GeV, top, and M^^ = 500 GeV, bottom) and background distributions 
after the Selection #1, #2 and #3 cuts. (Here, ^ = 100 fb^^.) 

#1-4 cuts, are given in Tab. IIV[ This summary clearly confirms the feasibility of the 
extraction of both signals after even less than 100 fb~^ of accumulated luminosity. 
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M^, = 200 GeV 



Cuts 


Ev. Signal 


Eff. % 


Ev. WZjj 


EfF. % 


Ev. tt 


EfF. % 


Ev. tilu 


EfF. % 


s/Vb 


1 


68.043(15) 


100 


5875.02(24) 


100 


99.699(3) 


100 


89.14(16) 


100 


0.87 


2 


68.043(15) 


100 


498.83(2) 


8.5 


5.3822(2) 


5.4 


19.38(3) 


21.8 


2.97 


3 


58.842(13) 


86.5 


10.5755(4) 


12.7 





0.8 


0.0667(1) 


2.2 


18.0 


4 


56.038(12) 


94.1 


4.4881(2) 


67.6 





56.4 


0.03047(5) 


64.8 


26.3 


M^^ = 500 GeV 


Cuts 


Ev. Signal 


Eff. % 


Ev. WZjj 


Eff. % 


Ev. ti 


EfF. % 


Ev. till' 


EfF. % 


S/^/B 


1 


73.668(31) 


100 


5875.02(24) 


100 


99.699(3) 


100 


89.14(16) 


100 


0.95 


2 


73.668(31) 


100 


498.83(2) 


8.5 


5.3822(2) 


5.4 


19.38(3) 


21.8 


3.22 


3 


68.833(29) 


93.4 


10.5755(4) 


12.7 





0.8 


0.0667(1) 


2.2 


21.1 


4 


46.337(20) 


66.0 


2.87857(1) 


7.1 





8.7 


0.00952(2) 


10.1 


27.6 



TABLE IV: Signal (M,,^ = 200 GcV at the top and = 500 GcV at the bottom) and background 
events per = 100 fb~^ and efficiencies Following the sequential application oF Selection #1-4 
cuts. 

IV. CONCLUSIONS 

We have analyzed the LHC discovery potential in the Z' and heavy neutrino sector of 
a (broken) U{1)b-l enlarged SM also encompassing three heavy (Majorana) neutrinos and 
found that novel signals can be established. The most interesting new signature involves 
three leptons (electron and/or muons), two jets plus missing transverse momentum coming 
from a Z' decay chain into heavy neutrinos. Various mass distributions (both invariant and 
transverse) can be used to not only extract the signal after a few years of LHC running, but 
also to measure the Z' and heavy neutrino masses involved. This is possible through DY 
production and decay via qq —>■ Z' —>■ Uh^h- In fact, for a large portion of the parameter space 
of our B — L model, the heavy neutrinos are rather long-lived particles, so that they produce 
displaced vertices in the LHC detectors, that can be distinguished from those induced by 
6-quarks. In addition, from the simultaneous measurement of both the heavy neutrino mass 
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and decay length one can estimate the absolute mass of the parent light neutrino, for which 
at present, only limits exist. 

This work has used a MC simulation based on a CalcHEP implementation of the B — 
L model. The analysis has been done at the parton level though we have verified that 
our results are stable against the implementation of typical ATLAS/CMS hadronic and 
electromagnetic calorimeter energy resolution effects. As benchmark scenarios of the B — L 
model we have chosen two that should be accessible at the LHC, having a Z' mass and 
fermion couplings not far beyond the ultimate reach of Tevatron and LEP and displaying 
two extreme relative conditions between the Z' and heavy neutrinos, that is, one with the 
latter produced at rest and the other highly boosted in the Z' direction. 
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APPENDIX: FEYNMAN RULES FOR HEAVY NEUTRINO INTERACTIONS 

In this Appendix, we list the Feynman rules involving the heavy neutrino of the B — L 
model considered. The intervening quantities are defined in the main text. 





W ^^Z 

V2e . e . ^ 

smau sm zau 



Hi 



H2 



— (—V2xy^C(x cos 2av + Mjy^ sin 2avSoc\ — (—V^xy^Sa cos 2ai, — M^^^ sin 2aiy 




Zb-l 




Zb-l 



g'l sin 2oLv 



g'l cos 2av 



where y — — , sm 2a,, — —2 



, cos 2av 



M, 
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